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ABSTRACT 

Accurate measurement of galaxy cluster masses is an essential component not only in studies of 
cluster physics, but also for probes of cosmology. However, different mass measurement techniques 
frequently yield discrepant results. The SDSS MaxBCG catalog's mass-richness relation has previously 
been constrained using weak lensing shear, Sunyaev-Zeldovich (SZ), and X-ray measurements. The 
mass normalization of the clusters as measured by weak lensing shear is > 25% higher than that 
measured using SZ and X-ray methods, a difference much larger than the stated measurement errors 
in the analyses. We constrain the mass-richness relation of the MaxBCG galaxy cluster catalog 
by measuring the gravitational lensing magnification of type I quasars in the background of the 
clusters. The magnification is determined using the quasars' variability and the correlation between 
quasars' variability amplitude and intrinsic luminosity. The mass-richness relation determined through 
magnification is in agreement with that measured using shear, confirming that the lensing strength 
of the clusters implies a high mass normalization, and that the discrepancy with other methods is 
not due to a shear-related systematic measurement error. We study the dependence of the measured 
mass normalization on the cluster halo orientation. As expected, line-of-sight clusters yield a higher 
normalization; however, this minority of haloes does not significantly bias the average mass-richness 
relation of the catalog. 

Subject headings: gravitational lensing:weak - galaxies: active ~ quasars: general - galaxies: clusters - 
methods:data analysis 



1. INTRODUCTION 

Galaxy clusters are powerful probes of cosmology, as 
they are the most massive collapsed structures in the 
universe. The mass function of galaxy clusters is sensi- 
tive to the matter density VLm, the amplitude of matter 
fluctuations us , and the density and evolution of dark en- 
ergy. Large catalogs of homogeneously selected clusters 
can therefore be used to constrain t h ese parameters (e.g . , 
IRozo et all IMot IWen et al.l [20T0I: ISehgal et al.l mm . 
Individual clusters can also provide important cosmolog- 
ical information. Hierarchical structure formation pre- 
cludes the formation of very massive clusters at high red- 
shift; it is unclear if recent measurements of distant large 
clusters are compatible with hie rarchical mass assembly 
in ACDM cosmo l ogy (see, e.g., IMortonson et al. 120111 : 
iFolev et al.l I20TI lYarvura et all 120111: IJee et al.ll201lD . 
Accurate cluster mass measurements are critical to each 
of these analyses. 

The mass of galaxy clusters can be estimated us- 
ing several methods. Weak lensing shear techniques 
involve measuring the total projected mass of a clus- 
ter by quantifying the gravitational distortion of the 
shapes of backgro und galaxies (for review s of w eak 
lensing the ory, see lBartelmann fc Schneiderl 120011 and 
ISchneider^ t al. 20Q1. For descriptions of s ome common 
weak lensing measurement techniques, see IKaiser et al.l 
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[iQQalRefregier fc Baconl[200a IMiUer et al.ll2007n. Ther- 
mal Sunyaev-Zeldovich (SZ: ISunvaev fc ZeldovichlflQTl 
analyses measure the change in the cosmic microwave 
background (CMB) due to inverse Compton scattering 
between electrons in the intra-cluster medium (ICM) and 
CMB photons. The amplitude of the SZ signal is de- 
termined by the electron number density and tempera- 
ture. By assuming a pressure profile for the clusters, one 
can convert the SZ signal into a cluster mass measure- 
ment (see Birkinshaw 1999, for a review). X-ray observa- 
tions also constrain the mass of a cluster, by measuring 
the thermal bremsstrahlung emission of electrons in the 
ICM; given assumptions about the dynamical state of 
the cluster, the X-ray lu minosity can be co nverted into a 
mass measurement (e.g.. lPratt et al.l[2009| ). These tech- 
niques often strug gle to agree when applied to the same 
clusters (e.g.. [Limousin et al.l l2010t IHuang et al.l 120111 : 
iMorandi et al.l 120111 ), alt hough some syste i ns do yield 
• good agreement (e.g..lZhang et al.ll2010l : llsrael et al.l 



very 

[25mlLCTchster et al.l ' l201llv r 

The MaxBCG catalog is an optically selected set of 
13,823 galax y clusters det ected in the Sloan Digital 
Sky Survey (jKoester et al.l I2007D . The average mass 
profile of the MaxBCG clusters has been extensively 
studied using weak l ensing shear (iSheldon et al.l 120091 : 
iJohnston etaH l2007t iMandelbaum et al.l l2008f K These 
studies have yielded precise measurements of the clus- 
ters' mass-richness relation, which connects the number 
of galaxies observed to be in the cluster (its richness) to 
the underlying halo mass. 

Recently, the SZ signal of the MaxBCG clus- 
ters has been measured by the Planck satellite 
(jPlanck Collaborationl[20"llal ) . The result implies a mass 
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normalization of the clusters that is > 25% lower than 
that determined using shear. This difference is signifi- 
cantly larger than the measurement errors quoted in the 
studies, and represents a puzzle with important implica- 
tions for both CMB and shear analysis techniques as well 
as for cluster physics and cosmology. 

The discrepancy in the mass-richness relations implied 
by the weak lensing shear and SZ analyses can be due 
either to systematic errors in the shear or CMB measure- 
ments, or to incorrect assumptions when comparing the 
results. The probability of errors in the basic measure- 
ments is small, due to the fact that both the shear and 
SZ results have been corroborated by additional investi- 
gations. 

The Planck results have been supported by CMB 
data from the Wil kinson Microwave A nisotropy Probe 
(WMAP) sateUite praper et al.ll201lD . Furthermore, 
an analysis of the SZ signal of clusters observed in X- 
ray data show good agreement between the mass esti- 
mates derived from the CM B and the X-ray analyses 
(jPlanck Collaboration! [2011 bO . It is important to note 
that while the X-ray and CMB data are entirely inde- 
pendent, they both measure the ICM gas and therefore 
the interpretations of the signals are subject to the same 
assumptions about cluster profiles and gas physics. 

The separat e wea k le nsing shear analy ses of 

[Johnston etldl ()2007[ ) and iMandelbaum et alj (|2008[ ) 
yield cluster masses that agre e at the lev e l of ^ 5%, as 
discussed in Appendix A of iRozo et "all ()2009f ) where 
the measurements are combined into a single relation 
with appropriate uncertainty, and small systematic 
differences between the two analyses are taken into 
account. However, the two works use similar shear 
measurement techniques (i.e. weighted second moments 
of the galaxies' flux distributions, corrected for shape 
respo nsivity using the results of iBernstein fc Jarvi^ 
120021) and similar background galaxy selections, making 
the analyses not completely independent. The lensing 
measurements are converted to cluster mass using the 
fact that the tangential shear at a distance R from the 
cluster center (on the plane of the sky) equals the critical 
density of the universe at the cluster's redshift, Ecrit, 
multiplied by AE: the average projected surface mass 
density inside a circle of radius R around the cluster, 
minus the average projected surface mass density at 
radius R. The shear measurements thereby yield cluster 
profile shapes, which are fit with Navarro -Frenk- White 
(NEW) radial profiles (|Navarro et al.l[T997[ ) to determine 
the cluster masses. (Only profile shapes, not ampli- 
tudes, are probed by weak lensing shear; magnification 
measurements do not suffer from this complication, 
which is known as the rnass-she et de generacy.) The fact 
that t he [Johnston et al.l ()2007| ) and IMandelbaum et aD 
(|2008f l results are consistent is an important check 
of the data analysis; however, because the analysis 
techniques are similar, the results are subject to similar 
systematic errors (for example, biases in galaxy shape 
measurement due to intrinsic prope rties of the galaxies 
and the image q uality, (see, e.g., iMassey et al.l 120071 : 
IBridle et al.l[20Tnfl V 

This work constrains the mass-richness relation of the 
MaxBCG cluste rs by measur i ng th e lensing magnifica- 
tion of quasars. ' Bauer et al.l ()201lD demonstrated that 
the lensing magnification of type I quasars can be mea- 



sured using the relation between the objects' variability 
amplitudes and their luminosities. A measurement of the 
variability of a quasar yields an estimate of its intrinsic 
luminosity, which can be compared to its observed lumi- 
nosity to provide a measurement of its magnification. 

Lensing magnification, like lensing shear, measures the 
total projected mass of the galaxy cluster acting as a 
lens. These two measurement techniques, however, are 
subject to very different systematic errors. This work 
is therefore an important test of the weak lensing shear 
analyses, checking that the procedures used to measure 
and interpret shape distortions do not suffer from a large 
systematic error that could dominate the discrepancy in 
the mass normalization of the MaxBCG clusters. Such 
an independent measurement of the cluster masses is a 
critical step in understanding and resolving the tension 
between the gas-based and lensing-based analyses. 

2. DATA 

We measure quasar variability using t he Palomar- 
QUES T RG-610 data set, as described in I Bauer et al.l 
(|2011f ) (hereafter called Paper I). The variability data 
are measured and calibrated identically to the data set 
in that work. In short, the Palomar-QUEST Variability 
Survey has observed 30,000 square degrees of sky, imaged 
multiple times in a broad, red, optical filter (RG-610). 
The survey reaches a depth of mag 19.5 in each exposure 
(under good conditions). After photometric calibration 
and strict quality cuts, the data exhibit a systematic er- 
ror of typically 2%. 

In the current work, the quasar sample is taken 
from the c atalog of quasar properties presented in 
iShen et al.l (l2011l i. The catalog con sists of 105,783 
quasa rs from SDSS data release 7 (jAbazajian et al.l 
l2009f l. and includes black hole mass estimates, lumi- 
nosity measurements at 1350A, 3000A, and 5100A, as 
well as flags marking, for example, broad absorption 
line (BAL) quasars and quasars observed in the ra- 
dio band. The variability-luminosity relation is seen 
in radio-quiet quasars, whose variability is thought to 
relat e to the behavior of accretion disk instabilities 
(e.g., 'CoUier & Peterson"2001; 'Vanden Berk et al."200i 
De Vrics et al. 2005; Rcngstorf et al. 2006; Wilhitc ct al] 
2008HBauer et al.ll2009bHMeusinger et al.il2011il . There- 
fore, we do not use in our analysis quasars that have 
observed radio emission, which implies the existence of a 
jet whic h can be an additi onal source of optical variabil- 
ity (e.g.. lKellv et al.|l2009[ ). We also eliminate from the 
sample those objects flagged as BAL quasars, as the out- 
flows from these objects may affect the observed variabil- 
ity properties. We investigate the effect of these cuts in 
section l4?2l We use the catalog's luminosities at 3000A 
as the luminosity measure in the variability-luminosity 
relation, as that wavelength exhibits continuum luminos- 
ity from the accretion disk (similar to the luminosity at 
25OOA, which was the measure adopted in Paper I). We 
use the mass estimates from the catalog as well, when 
normalizing the variability data as described below. 

3. METHOD 

We use the technique presented in Paper I to mea- 
sure the lensing magniflcation of type I quasars using 
the quasar variability-luminosity relation. We briefly de- 
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scribe the method below. For further explanation, 
refer the reader to Paper I. 



we 



3.1. Measuring the Variability- Luminosity Relation 

The basis of the magnification measurement is the em- 
pirical correlation between variability amplitude and lu- 
minosity that has been seen in large ensembles of type I 
quasars (jVanden Berk et al.l [200l iBauer et all l20'09b[ ). 
This relation allows a measurement of the variability 
amplitude of a quasar to act as a measurement of the 
quasar's intrinsic luminosity. The measured luminosity 
of the object, divided by this intrinsic luminosity, yields 
an estimate of the quasar's magnification. Each mag- 
nification measurement has a large error, corresponding 
to the large scatter in the variability-luminosity relation. 
We can average the magnification measurements of many 
quasars, however, to obtain a statistically significant re- 
sult. 

We measure the variability of quasars using the statis- 
tic 

V = v'(Am)2 - a2 (1) 

which is similar to the structure function. Am is the 
magnitude difference between two measurements of a sin- 
gle quasar, and a is the error on those photometry mea- 
surements. V is measured for each pair of magnitude 
measurements of a quasar; N observations of one object 
will yield A''(A'' — 1) measurements of V. 

Quasar variability amplitude has been seen to depend 
on a number of factors besides luminosity, for example 
the time lag between measurements, the quasar's esti- 
mated black hole mass, and wavelength of the observa- 
tions ('e.g..lVMi dcn Ber k et al.ll2004 lDe Vries et al.ll2Q05l: 



Wilhite et al.' 2005; Wilhite et al.' 2008; Bau er et al.l 
2P09b; MacLeod et al. 2010; Mcusingcr et al. 20111). In 
order to measure a tight correlation between the vari- 
ability and luminosity, it is necessary to take into ac- 
count the variability amplitude's dependence on these 
other parameters. This is done by dividing the quasars' 

V measurements into bins according to black hole mass 
M, wavelength of observation A, time lag t, and lumi- 
nosity L. Each multi-dimensional bin yields an average 

V value; the difference in these average Vs is due to the 
dependence of V on mass, wavelength, time lag, and lu- 
minosity. For each choice of (M, A, r) bins there exists a 

V — L relation whose amplitude is related to the depen- 
dence of V on M, A, and r, but whose slope is due to the 
relation between V and L. We determine a zero point 
for each (M, A, t) bin, additive in log(y), such that the 
different bins' Vs are normalized to the same amplitude. 
The variability's dependence on M, A, and r is there- 
fore removed, and the normalized data demonstrate how 

V scales with L. Explicitly, the variability- luminosity 
relation is of the form 

log(y) - Co + SC{M, A, r) - a X log(L„,eas). (2) 

where Co is a constant to be fit using the data, and 
5C{M, A, r) are the zero points that normalize the dif- 
ferent bins' Vs to the same amplitude. We then define 
the normalized variability Vnorm such that 



log(K, 



log{V) ~ SC{M, A, r) 



(3) 
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Table 1 

Bin limits used in determining the normalization constants. 
Measurements with values outside the limits are not considered. 
Units of time lag t: days; black hole mass M: 10* X Mq, 
luminosity L: 
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observe a quasar at redshift z is simply Ar,g610 x (1 + 
where Arg610 is the central wavelength of the RG610 
filter. Each measurement includes flux from a range of 
rest frame wavelengths because the filter is very broad. 
However, the observed wavelength range scales uniquely 
and monotonically with redshift. Therefore, binning the 
quasars in redshift z is equivalent to binning them in 
A. If quasar variability evolves with redshift in addi- 
tion to depending on observed wavelength, then using 
single-filter data and binning in either z or A confiates 
the two effects. Recent studies measure no significant de- 
pendence of variability on redshift ( MacLeod et al. 201(1 
iMeusinger et al.ll2011f) , implying that by binning in z we 
are primarily capturing the change of variability ampli- 
tude with rest frame wavelength. 

The bin limits in luminosity, black hole mass, redshift 
(in place of wavelength) , and time lag are given in Table 

m 

3.2. Data Cuts 

We make similar quality cuts to the data as in Pa- 
per I. The calibrated Palomar-QUEST Survey data in- 
clude 819,679 measure ments of 13,782 ob jects from the 
DR7 quasar catalog of iShen et al.l ()2011|) . after remov- 
ing BAL quasars and objects with radio detections. We 
further wish to eliminate from the sample quasars that 
show unusual variability behavior, as they are less likely 
to follow the variability-luminosity relation seen for the 
ensemble. For example, a quasar may have a jet which 
is not bright enough to be noted in the catalog, but still 
contributes to the optical variability. Because the vari- 
ability mechanism for this object would differ from that 
of the ensemble, we would not expect it to follow the 
same variability-luminosity relation. We therefore elimi- 
nate objects that exhibit variability that is significantly 
larger than typically seen in type I quasars. In particular, 
for each object we calculate the structure function 



SF ^ y/< {Amy 



> - < o-^ > 



(4) 



Because the Palomar-QUEST data use a single optical 
filter, the central rest-frame wavelength A at which we 



where angular brackets indicate the mean. Am and cr 
are as in equation [TJ The structure function SF is cal- 
culated twice for each object: once including all mea- 
surement pairs Am with time lag less than 10 days in 
the rest frame, and once including all measurement pairs 
Am with rest frame time lag greater than 100 days. If the 
short timescale SF is larger than 0.1 or the long timescale 
SF is larger than 0.5, the object is discarded from the 
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sample. These cut values are motivated by the variability 
of lar^e ensemb l es of ty pe I quasars and blazars measured 
in iBauer et"all (|2009a[ ). We thereby eliminate 1237 ob- 
jects, leaving 12,545 quasars with 737,730 measurements. 
Note that this cut also eliminates objects that may be 
poorly measured or calibrated in the survey. 

Many objects are cut because they do not lie inside the 
bin limits given in Table [TJ For example, there exists a 
long, low-level tail of measurements with time lags larger 
than 400 days. Such tails are difficult to calibrate well, 
and should be eliminated from the sample. After binning 
the data in time lag, redshift, mass, and luminosity we 
are left with 463,437 measurements of 8499 quasars. 

To further eliminate poorly measured data from the 
sample, we cut any measurement pairs Am that exhibit 
variability greater than one magnitude, as this is unex- 
pected in type I quasa rs over time lags studied in this 
work fsee IBauer et al.l r2009a.' ) . This cuts only 17 mea- 
surement pairs. 

As in Paper I, we consider only measurement pairs for 
which we observe significant variability, such that V > a 
in equation [TJ This cut eliminates much of the data, 
but is necessary in order to ensure that the variability 
measurements are not dominated by noise. The result- 
ing data set includes 166,101 measurements from 7713 
quasars. 

When calculating the zero points for each multi- 
dimensional bin of mass, wavelength, and time lag, we 
insist that there exist 30 measurements in each bin such 
that the zero point is well-determined. If a bin has fewer 
members, it is not considered in the analysis. This crite- 
ria leaves 165,845 measurements of 7712 quasars. 

Finally, because not all quasars have the same num- 
ber of measurements, and because a magnification mea- 
surement is calculated from each measurement pair (see 
equation|6]), a minority of quasars with a large number of 
measurements may have a disproportionate influence on 
the results. We therefore include in the analysis no more 
than 20 measurement pairs from each quasar. This cut 
produces our final data set, with 97,247 measurements 
of 7712 quasars. 

The effe cts o f these cuts on the results are investigated 
in section 

3.3. Measuring Quasar Magnification 

After the V measurements are normalized, the 
variability-luminosity relation is measured by fitting the 
parameters a and C in the equation 

log(Korm) = Co - a X log(Lincas)- (5) 

where imeas is the quasar luminosity at 3000A, and the 
relation is fit using the entire quasar sample. Then, as- 
suming that the mean magnification of the entire quasar 
sample is unity, the magnification is given by 



in 



IQCo 



1/q 



(6) 



Equation [S] is used to calculate a magnification estimate 
from each V measurement in the data set. We can as- 
sume that the mean magnification of our quasar sample 
is indeed unity, as it has been shown in ray-tracing sim- 
ulations that the mean magnification across the sky is 
unity (,Takahashi et al.,.2011, ). and the current data set 



spans thousands of square degrees of sky and therefore 
should cover a fair sample of large scale structure. 

The error in a magnification measurement reflects the 
scatter in the variability-luminosity relation, which is 
due to measurement uncertainties and also to details of 
quasar variability that are not well understood. The er- 
ror on each measurement is larger than the magniflca- 
tion signal; therefore, we must average many magniflca- 
tion measurements in order to obtain a significant result. 
In order to combine measurements with similar expected 
magnification, we bin the measurements by their distance 
(scaled by the cluster radius) from the nearest galaxy 
cluster in the MaxBCG catalog. 

3.4. Constraining Galaxy Cluster Profiles 

The MaxBCG catalog contains 13,823 clusters de- 
tected in the SDSS, with redshifts between 0.1 and 0.3 
and masses greater than roughly 10^'* Mq. The catalog 
is estimated to be 90% pure and 85% complete in its 
redshift range and 7500 square degree footprint. 

To measure the average profile of the clusters in the 
catalog, we average the magnification signal of quasars 
close to the clusters. To rescale the clusters' magnifica- 
tion signal so that we can effectively stack clusters with 
very different richness, we consider the magnification as 
a function of distance from the line of sight to the cluster, 
divided by i?200 of the cluster, where i?2oo is the radius 
inside which the density is 200 xpcrit, the critical density 
of the universe at the redshift of the cluster. In order to 
determine i?200i ^^nd to calculate the magnification ex- 
pected from the clusters, we must assume a cluster mass 
model. We assume t hat each cluster has a NFW profile 
(|Navarro et al.lll997[) : 



p{r) 



^cPcrit 



ir/n)[l + ir/rs)Y^ 



(7) 



where c is the concentration parameter of the profile, 

Sc = ^ in(i+c)-c/(i+c) ' ^^^'^ = ^200 /c. The total 
mass inside i?200 is Af2oo- 

In order to specify a unique model for the cluster 
masses, we must adopt a mass-concentration relation and 
a mass-richness relation. Given these assumptions, the 
richness parameter of each cluster in the MaxBCG cata- 
log fully determines the mass profile model. We assume 
the mass-concentration r elation from lDuffv et al.l ([2008): 



c(M2oo,z) = 5.71 



200 



-0.084 



(1 + ^)- 



(8) 

The NFW magnification profile can be calculated an- 
alytically using the p rofile parameters, as explained in 
ISchneider et"al1 (l200l . 

The mass richness relation is the quantity we wish to 
constrain. We a dopt the form of the relation used in 
iRozo et all (|2iD09t ): 

M20o(iV20o) = i^^y lO'^h-'M^ (9) 

where A'^200 is the number of red sequence galaxies ob- 
served to be in the cluster that are above a limiting 
brightness and within a certain radius, and is provided 
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Figure 1. Magnification versus transverse distance to tlie near- 
est cluster (scaled as R/R20o)- Stacked signal using all MaxBCG 
clusters. The solid line shows the expected magnification profile, 
calculated assuming a mass-richness relation of the form given in 
equation |9] with A = 0.49 and B = 0.95, which has a 16% proba- 
bility of being consistent with the data. The bottom panel shows a 
profile measured after randomizing the measurements, as described 
in the text, which is consistent with zero. 

in the MaxBCG catalog fsee lKoester et alj [20071 for de- 
tails). 

Our fitting procedure involves choosing A and B, cal- 
culating the expected magnification of the quasar sample 
in logarithmic bins versus radius, and comparing it to the 
data in order to calculate the probability that the data 
are consistent with the model. We renormalize the ex- 
pected magnification values such that the mean over the 
entire quasar sample is equal to one; this typically shifts 
the values by only about 0.4%. We choose A and B from 
a grid of points with spacing 0.01 in each dimension, from 
{A,B) of (0.1,0.8) to (0.8,1.3). 

4. RESULTS 

An example plot of the expected and measured magni- 
fication, versus scaled distance from the nearest cluster, 
is shown in the main panel of Figure [TJ The bins' ex- 
pected magnifications are connected by lines for clarity. 
The values of A and B are 0.49 and 0.95, respectively. To 
avoid the fit quality from being dominated by properties 
outside of the main cluster halo we consider five bins, 
out to radius 3 x -R200- This particular combination of A 
and B provides the best fit to the data, with a chi square 
of 8.0 with 5 degrees of freedom, corresponding to 16% 
probability of being consistent with the data. 

The sources and detailed properties of the scatter in 
the variability-luminosity relation, which is the main 
source of error in the magnification measurements, are 
not well understood. Therefore, we calculate the magni- 
fication errors empirically using the distribution of mea- 
surements. We calculate the error for each radial bin us- 
ing bootstrap resampling, with 500 subsamples per bin. 
Specifically, for each bin we randomly choose N objects 
from the data in that bin, where N is the number of data 
points in the bin. Because we choose randomly with 
replacement, the bootstrap sample is different from the 
original data sample in that some values are missing and 
others are repeated. We construct 500 such bootstrap 
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Figure 2. Mass-richness relation parameters A versus B, with 
gray scale corresponding to the probability of agreement with the 
data. The contours mark 68.3%, 95.5%, and 99.7% confidence re- 
gions around the best -fit values. Boxes repre sent (from t o p to bot- 
tom) th e 1(T regions of lRozo et al.l l|2009l ) and lPlanck Collaboration! 
pOllal ). 

samples, and take the error to be the range that encloses 
68% of the samples. We have tested that the errors are 
insensitive to the exact number of bootstrap subsamples, 
and that randomly shuffling the magnification measure- 
ments (i.e. replacing each measurement with another 
from the data set in a random manner) produces a profile 
consistent with zero to within la. An example instance 
of the profile using randomized measurements is shown 
in the bottom panel of Figure [1] 

4.1. Mass- Richness Relation 

Figure [2] shows the mass-richness parameters A ver- 
sus B, with the gray scale corresponding to the proba- 
bility that the mass model is consistent with the data. 
The contours mark the 68.3%, 95.5%, and 99.7% (in 
order of decreasing line thickness) confidence regions 
around the best-fit parameters. The upper and lower 
boxes delineate the la result s of iRozo et al.l ()2009[ ) and 
IPlanck Collaboration! (j2011af) . respectively. The best- fit 
model parameters, with marginalized 68% statistical er- 
rors, are {A, B) = {OAdto fi^O-Q^to fi)- The best-fit 
values have a 16% probability of agreeing with the data, 
and yield the profile shown in Figure [TJ 

Rozo et al. (2009) combine previous shear analyses to 
conclude {A, B) = (0.48 ± 0.07(stat) ± 0.06(sys), 1.13 ± 
0.09(stat) ± 0.()5(sys)). These systematic errors include 
the effects of cluster miscentering on the average shape 
of the haloes, due to incorrect identifica tion of the cen- 
tral g alaxies as discussed in detail in [Johnston et al.l 
(pool . The Planck SZ measurements are translated 
to {A, B) values using the fit results from Table 2 of 
IPlanck Collaboration! ([2011aD (using the lRozo et alll2009l 
mass calibration) for the Isoo — ^^200 relation, where Isoo 
is the SZ signal: the integral of the electron pressure 
over the volume of the cluster, times a geometric factor, 
scaled to a redshift of zero so that the measurements of all 
of the clusters are directly comparable. The conversion 
between i^oo and Msp n, as put forth in 'Arnaud et al.l 
(i2010.) and adopted in [Planck Collaboration. (,2011al) is 
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such that Ysoo oc AfjQQ . M^qq is the mass of the cluster 
inside a radius within which the average matter density 
is 500 times the critical density of the universe at the 
redshift of the cluster. Making the simplifying assump- 
tion that the conversion from M500 to M200 does not in- 
troduce error into the mass-richness relation parameters, 
the fractional error on the M500 — N200 parameters from 
the SZ analysis will be the same as the fractional error 
on the equivalent M200 — -^200 relation parameters. We 
can therefore propagate the errors from the Y^qq — N2Q0 
relati on to approximate err o rs in a n M200 — -^200 rela- 
tion. iPlanck Collaboration! ()2011af) state that their re- 
sults ar e in agreement with the mass-richness relation 
given in lJohnston et al.l (|2007f ) , except wi th an amplitude 
smalle r by 25%. We therefore use the Joh nston et al.l 
(j2007D results to determine the Planck-based value of 
A = 0.19 ± 0.03. The Planck-based B, calculated from 
their parameter a, is 1.22 ± 0.04. 

4.2. Ejfects of Data Cuts 

To test that the cuts made to the data set do not bias 
the cluster mass measurements, we repeat the analysis 
using data samples with modified selection criteria. 

In the main analysis, we eliminate from our sample 
quasars detected in the radio band by the FIRST Survey 
(jBecker et al.lll995l ). which are flagged in the value-added 
quasar catalog. Because jets produce radio emission and 
also affect quasar optical variability, we expect those 
quasars with radio emission to introduce extra scatter 
into the variability-luminosity relation, creating noise in 
the cluster mass measurements. Because there are bound 
to be quasars in the sample with radio fluxes just be- 
low the FIRST detection limit, it is important to ensure 
that radio emission is not associated with strong biases 
in the mass measurements. To explicitly test the effects 
of radio-detected quasars on the mass measurements, 
we repeat the above analysis while including quasars 
with radio fluxes measure d to be l ess th an 10 mJy at 
6 cm (as given in the Sh en et all ()2011[ ) quasar cata- 
log). This comprises 52% of the radio-detected quasars 
in the catalog. We find that similar results when includ- 
ing the radio detections. Specifically, with 68% confi- 
dence, iA,B) = (0.50t|]:[!^,0.98l:^:[5^). The inclusion of 
the radio-detected quasars enlarges the confidence region 
to higher mass normalizations, as expected by the fact 
that radio-loud AGN show higher variability than those 
that are radio-quiet. The effect, however, does not sig- 
nificantly bias the results. 

Including only variability measurements V that are 
larger than 2a rather than la, eliminating highly vari- 
able objects with long- and short-term SF 25% larger or 
smaller than the standard cut, and allowing only 10 mea- 
surements per objects rather than 20, each changes the 
best-fit {A, B) by less than the 68% errors. Eliminating 
the typically least-variable measurements, by discarding 
the first two bins in time lag, similarly has no significant 
effect on the results. 

Including the 5258 AGN marked as BAL quasars in 
the value-added catalog decreases the best-fit A and B 
along with the probability of the best-fit model repre- 
senting the data. The best-fit region lies along the same 
degeneracy direction as is visible in Figures [2] and [3l al- 
though the highest probability values are at the low-i? 



range of the parameter space. Because including known 
BAL quasars in the sample significantly lowers the best- 
fit [A, B) values, unidentified BAL quasars remaining in 
the sample may be artificially lowering the best-fit mass- 
richness parameters for our main quasar sample. Cor- 
recting for such conta mination would ser ve to improve 
the agreement with the iRozo et"al] (|2009f) mass-richness 
relation. 

To examine the sensitivity of our results to the chosen 
concentration relation, we repeat the analysis using an 
alternate relation, from Mandclbaum et al. (2008); 

c[M2,,,z) = ^( ^^^2^-^ (10) 

This relation changes the results by less than the 68% 
errors. 

4.3. Effects of Cluster Miscentering 

The MaxBCG catalog is known to suffer from miscen- 
tering uncertainties. The cluster centers are given as the 
location of the brightest cluster galaxy (BCG). In some 
cases, the BCG may be offset from the true halo cen- 
ter. More often, the BCG is misidentified. Such errors 
in the determination of the halo center will lead to er- 
rors in the stacked halo profile, in particular making the 
stacked profile shallower than expected. This effect can 
show up as an artificial change in co ncentration depend- 
ing o n the innermost radius probed ([Mandelbaum et al.l 
HOQD- [Johnston et al. (200% studied this problem in 
detail using simulations, resulting in a mass-dependent 
probability that a cluster is offset from its true center by 
a given amount. 

The fact that our results are robust to changes in the 
cluster concentration relation implies that they may not 
be strongly affected by miscentering. To investigate this 
further, we implemeii t the miscentering prescription of 
[Johnston et al.l (|2007[ ). offsetting the centers of a ran- 
dom minority of the clusters. This procedure likely over- 
estimates the true miscen tering properties of the catalog 
(jMandelbaum et "an '2008') . It is not meant to correct for 
the cluster center misidentification, since the randomly- 
chosen offsets will almost always bring the considered 
coordinates farther away from the true cluster center. 
Instead, applying this prescription is meant to show the 
effects of such center shifts on the results. The proce- 
dure adds noise to the results, while moving the best-fit 
region along the direction of (A, B) degeneracy towards 
lower best-fit values. The precise best-fit region is sen- 
sitive to the random choice of cluster offsets; different 
realizations of the procedure yield different results. Be- 
cause the best-fit values consistently lie along the line of 
{A, B) degeneracy, miscentering is un able to bring our re- 
sults in line with the predictions of Planck Collaboration! 
(l2011aD . If it were possible to correct for cluster miscen- 
tering, our best-fit A and B would both increase. Such 
behavior would not ease the tension between the weak 
lensing and SZ results. 

4.4. Effects of Cluster Orientation 

Brightest cluster galaxies (BCGs) have been observed 
to have major axes that are typi cally aligned with that 
of their host dark matter haloes ([Hashimoto et al]|2008l : 
iNiederste-Ostholt et al...2010t ). The morphologies of the 
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Figure 3. Mass-richness relation parameters A versus B, with 
gray scale corresponding to the probability of agreement with the 
data, for the "face-on" cluster sample. The contours mark 68.3%, 
95.5%, and 99.7% confidence regions around the bes t-fit values. 
Boxes represent (from top to bott om) the la regions of lRozo et al.l 
(|2P09) and Planck Collaboration pOllal ). 



Figure 4. Mass-richness relation parameters A versus B, with 
gray scale corresponding to the probability of agreement with the 
data, for the "line-of-sight" cluster sample. The contours mark 
68.3%, 95.5%, and 99.7% confidence regions around the best-fit 
values. Boxes re prese nt (from top to bottom) the la regions of 
IRozo et al.1 l(2009l ) and lPlanck Collaborationl l|2011al ). 



clusters, as opposed to their axes of orientation, are not 
seen to be related to the morphologies or orientations 
of the BCGs. Because of this connection, clusters with 
BCGs that appear elliptical on the plane of the sky are 
more likely to be oriented face-on than those with BCGs 
that appear spherical (because an apparently-spherical 
BCG may be an elliptical galaxy with major axis close 
to t he line of sight). 

Marrone et al.l (|2011l ). in comparing the weak lensing 
and SZ signals from 18 galaxy clusters, measured a scat- 
ter around the mean relationship between shear-based 
mass and SZ signal that was dependent on cluster mor- 
phology and orientation. The clusters' morphology was 
determined using X-ray data, and undisturbed clusters 
yielded systematically larger weak lensing masses than 
disturbed clusters with similar SZ signal. The clusters' 
orientation was approximated by the ellipticity of their 
BCGs; those with axis ratio larger than 0.85 (i.e., clusters 
preferentially oriented along the line of sight) also had 
systematically larger weak lensing masses for a given SZ 
signal. These results agree with the expected relation- 
ship between halo orientation and bias in cluster weak 
lensing mass measurements (Corlcss & King 2007). 

Because X-ray measurements are not available for the 
majority of the MaxBCG catalog, it is difficult to esti- 
mate the morphology of our cluster sample. The elliptic- 
ities of the BCGs are available from the SDSS CasJobs 
serveiQ; we use these as a proxy for the orientation of 
the MaxBCG clusters. We repeat our analysis using only 
those MaxBCG clusters for which the BCG has axis ra- 
tio less than 0.85. By matching the MaxBCG cluster 
central position with the SDSS DR7 galaxy table, with 
a spatial matching tolerance of 1" , we find 10,248 BCG 
galaxies with axis ratio, calculated by fitting an elliptical 
deVaucouleurs profile to the r band SDSS data, of less 
than 0.85. The mass-richness parameters determined for 
this "face-on" sample is shown in Figure |3l The results 
agree within la with those of the whole sample, although 



^bcst ^bcst Pbcst Prozo Ppianck 

Total Sample 0.49+q;^jJ; 0.95+^;^^ 0.16 0.11 7.0e-06 

Face-on 0-i'2tom+ 0.81^0.0?+ 0.14 0.08 1.5e-05 

Line-of-Sight 0-69to;o3_). 0.83to;o3+ 0.64 0.04 5.7e-08 



Table 2 

Best-fit mass-richness parameters A and B and the probability 
that the magnification model with these parameters represents 
the data. Prozo and Ppianck a^re the t otal probab ility within the 
la error regions of the results from lHozo et ahl 820091 ) and 
Planck Collaboration (2011a), respectively, normalized to the 
probability in the region from (A, B) (0.1,0.8) to (0.8,1.3). The 
"face-on" and "line-of-sight" subsamples' error regions extend 
outside of^ our parameter space. 

the best-fit region extends to significantly lower A and 
B values. 

To compare with the face-on cluster sample, we con- 
struct a MaxBCG subsample including the 10,248 clus- 
ters with the roundest BCGs. There is substantial over- 
lap between this subsample and the "face-on" one used 
above; this subsample includes clusters with BCG axis 
ratio greater than ~ 0.68. Nevertheless, this "line-of- 
sight" subsample can serve as a comparison to highlight 
the effects of cluster orientation on the measurement. 
The mass-richness fit results for the "line-of-sight" sub- 
sample are shown in Figure ID The best-fit relation in- 
dicates significantly higher mass normalization for these 
"line-of-sight" clusters than for the whole sample and the 
"face-on" subset. Furthermore, the probability that the 
model represents the data is much higher for the best-fit 
parameters in the "line-of-sight" subsample (64%) than 
in the "face-on" subsample (14%), implying that corre- 
lated structure in the plane of the sky in face-on clusters 
decreases the quality of a spherical NFW fit. We note 
that the subsamples' best-fit regions extend beyond our 
explored parameter space, although the probed region is 
sufficiently large for the purpose of comparing the sub- 
sets' general behavior. 



^ http://cas.sdss.org/astrodr7/ 



5. DISCUSSION AND CONCLUSIONS 
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We constrain the mass-richness relation of the 
MaxBCG galaxy cluster catalog by measuring the lens- 
ing magnification of background quasars. The magnifica- 
tion is determined by measuring the quasars' variability 
and using the correlation between their variability and 
luminosity to constrain the objects' intrinsic luminos- 
ity. Comparing this intrinsic luminosity to the measured 
value yields a measurement of the quasar's magnification. 
Because each magnification measurement has low signal 
to noise, we stack the measurements to generate a single 
radial magnification profile as a function of distance from 
a cluster. 

Each galaxy cluster is assumed to have a NFW density 
profile, with a mass that is related to its richness iV2oo 
that is given in the MaxBCG catalog. The mass-richness 
relation is assumed to be of the form given in equation O 
Assuming values for the relation's normalization A and 
slope B, the magnification expected for the quasars, due 
to the galaxy clusters, can be calculated. The probabil- 
ity of the assumed mass-richness parameters is taken to 
be the agreement between the measured magnification 
and the expected values calculated using the clusters' 
mass- richness relation and the lensing properties of NFW 
haloes. These probabilities are used to determine the 
level of agreement between the data and the model, and 
bet ween the data and the result s of iRozo et al.l (|2OO90 
and iPlanck Collab oration! (|2011aD . The results are sum- 
marized in Table [21 

The best-fit and 68% statistical errors for the mass- 
richness parameters are {A,B) = (0.49lo;gi, O.95io;g5). 
Assuming that the model used to calculate the cluster 
magnification is accurate (i.e. that the clusters follow 
NFW profiles which obey a mass-richness relation param- 
eterized by A and B), we find ff% agreement between 
the magnification results and those of Rozo et al. ( 200!^. 
In contrast, we find 7e-04% agreement between the mag- 
nification data and the mass- richness relation im plied 
by the SZ results of IPlanck Collaboration! (|2011a! ). The 
magnification results thereby corroborate the weak lens- 
ing shear measurement of the MaxBCG mass-richness 
relation, and strongly disfavor the lower mass normal- 
ization implied by SZ and X-ray analyses. The data, 
however, are not perfectly represented by the model; the 
unnormalized best-fit probability that the data are rep- 
resented by the model is only 16% for the best-fit pa- 
rameters. We have checked the response of the results to 
cuts in the data, and to cluster miscentering. Systematic 
errors such as cluster miscentering and inclusion of het- 
erogeneous AGN in the sample (such as BAL quasars) 
serve to decrease A and B simultaneously, along a line of 
degeneracy visible in Figures [2] and [3l These effects are 
larger than the statistical errors; however, correcting for 
them would not bring the results closer to those of the 
SZ and X-ray analyses. 

The normalization of the mass-richness relation de- 
pends on the orientation of the galaxy clusters with re- 
spect to the line of sight. We split the MaxBCG cata- 
log into two overlapping sets, each of size ~ 75% of the 
whole catalog, based on the measured ellipticity of their 
BCGs which can be used as a diagnostic for the clus- 
ter halo orientation. The "face-on" subset (with BCG 
axis ratio less than 0.85) shows behavior similar to the 
whole sample, with consistent best-fit values although 
the 68% confidence region extends along a line of de- 



generacy to smaller A and B values. The "line-of-sight" 
subset (with BCG axis ratio greater than 0.68), which 
has identical statistics (and shares ~ 50% of its objects 
with the "face-on" subset) has a best-fit mass normaliza- 
tion A that is 40% higher than that of the whole sample. 
This is in line with predictions from numerical simula- 
tions for the lensing mass bias due to assuming spherical 
NFW profiles for prola te, line-of-sight oriented clusters 
(jFeroz fc Hobsonl!2011l) . Although the minority of clus- 
ters oriented close to the line of sight fit to a high mass 
normalization, they are not responsible for the disagree- 
ment between the main result and the SZ best-fit region, 
si nce both subsamples a r e incom patible with the results 
of !Planck Collaboration! ()2011a! ). 

The magnification-based measurement of the MaxBCG 
galaxy cluster mass-richness relation presented in this 
work is based on an analysis technique with very different 
systematic errors from those of weak lensing shear. Be- 
cause the mass-richness relation of the MaxBCG cluster 
catalog is consistent when measured using lensing shear 
and lensing magnification, the discrepancy between the 
lensing-based cluster masses and the SZ and X-ray ICM- 
based masses is unlikely to be related to a systematic 
problem with the shear measurement techniques. In- 
stead, the cause is likely a physical assumption about the 
galaxy clusters, made by one or both analyses. One as- 
sumption made by all of the analyses is that the clusters 
are well-represented by a single, spherical profile shape, 
while in fact the clusters are triaxial. Weak lensing shear 
measurements are expected to have biase s due to cluster 
orientation fsee !Becker fc Kravtsovll20Tol and references 
therein), moreover due to the combination of orienta- 
tion and halo shape (i . e. whether it is prolate or oblate) 
(jCorless fc Kind 120081: ll^roz fc HobsonI l20Tl[) . O ur re- 
sults indicate that although the orientation of the tri- 
axial haloes does affect the mass measurement, it does 
not significantly bias the average mass-richness relation 
of the MaxBCG catalog. Further investigation is needed 
into assumptions in the lensing and SZ analyses in order 
to resolve the significant discrepancy between the two 
methods' measurements of the MaxBCG cluster catalog. 
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